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ABSTRACT 


In the current study, it was tried to describe a method for determining 
thermal characteristics of integrated micro-circuits to identify thermal 
parameters of multidisciplmary (thermal-electric) behavioral models. 
The problem is solved on the example of high-frequency pulse voltage 
converters. A solution was proposed to refine the mimimum structure 


of the thermoelectric model based on an iterative least squares method using 

the Levenberg-Marquardt algorithm, as well as a graph of the spectral 
Keywords: den-sity of time constants. This made it possible to reduce the influence 
Integrated circuits of the filtering factor in the deconvolution Operation when building 

a thermal model using the structural function of the thermal characteristic 
Matlab transition. Also, the results obtained can be used to build integrated circuits 
P-n-junction (IC) behavioral models, taking into account the thermal processes 
Thermal resistance occurring in them. 
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1. INTRODUCTION 

Integrated circuits (IC) of high-frequency switching voltage converters (HFSVC) are widely used 
in power supplies of any electronic devices that make up the onboard systems of aircraft and spacecraft, 
household apphances, computers, mobile gadgets, etc. The area of their use as regulators of constant voltage 
is also constantly expanding due to the development of areas where linear constant-voltage regulators 
of the compensation type were used previously. 

Under modern conditions, shortening the development time of electronic devices is impossible 
without modeling. When developing the design of power electronics devices, one of the important tasks of 
modeling is to take into account the mutual influence of the thermal and electrical characteristics of 
semiconductor devices and microcircuits, which affects the reliability and performance of the electronic 
device as a whole. Therefore, the practice of applying complex multidisciplmary models is expanding. 
Usually, finite element models with distributed parameters [1-5] or with lumped parameters [6-9] are used. 
The first type of models allows modeling, taking into account many design features, but it requires a larger 
amount of input data for model compilation as well as a high level of scientific competence of the developer, 
which is not often observed in engineering practice. In this regard, the second type of the model has a number 
of advantages. Such models of integrated circuits can be successfully implemented in the popular electronic 
industry languages SPICE, VHDL-AMS, Verilog-A, etc. [9]. 

The full-fledged synthesis of multidisciplmary behavioral models of IC, based on the analysis of 
the electrical circuit and the design of the core of the chip, is an independent task and is beyond the scope of 
this article. The authors' task is to find and substantiate a method for determining the IC thermal 
characteristics, which allows it to be used in the synthesis of multidisciplinary behavioral models of high- 
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frequency impulse-voltage chip circuits, taking into account the features inherent in a particular HFSVC 
design. This paper solves the problem of obtaining a model describing the HFSVC -thermal properties with 
respect to the experimentally measured thermal transient response (TTR) the dependence of the crystal 
(junction) temperature of the circuit on time TJ (t). 


2. LITERATURE REVIEW 

To determine the TTR in the standard [10], it is recommended to periodically connect 
the semiconductor device under test to the nominal mode with a frequency of less than 66.7 Hz, to connect 
a small measuring current between pulses in the pauses of 100 us as well as measure the UF forward voltage 
at the p-n-junction. Using the UF (TJ) calibration, it is possible to obtain a set of TJk transition temperature 
values at the corresponding tk time points or a graph of the nonlinear TJk (tk) dependence. 

In the experimental TTR determination, direct measurement of the crystal temperature is impossible. 
The placement of the thermocouple in the instrument case (in close proximity to the crystal) or 
on the case [11] leads to a distortion of the real picture of the thermal processes taking place in 
the semiconductor device. In this regard, the methods of indirect determination of the crystal temperature 
of the device, based on the control of temperature-dependent parameters (TDP) inherent in the device under 
test, are actively used. The direct current voltage drop on the TJ p-n junction, the transconductance 
of the transistor S, reverse current (IRR) and the characteristics of the reverse pn junction recovery process 
(tRR, IRR,max), signal delay-propagation time tD [12-14] are usually used as IC TDP. When choosing 
a TDP HFPVC circuits, the following conclusions can be made. Due to the lack of connection 
from the external leads to the input circuit of the power transistor as part of the chip, it is not possible to use 
the transistorslope as the TDP. 

The use of tRR, IRR, max, tD dynamic parameters as a TDP makes it possible to experimentally 
determine the TTR of microcircuits directly in the nominal mode of the microcircuit; - the measurement of 
heating power requires the dynamic measurement of the microcircuit power consumption and the value of its 
“transit component” attributable to passive elements (choke, capacitors, external power diode) as wellas load 
included in a typical circuit switching on HFSVC. 

If TJ is used as a TDP, then usually difficulties in determming the TTR in the nominal mode of 
the chip are also associated with the implementation of heating power dynamic measurement, and 
additionally with the HFSVC tss start-up time, which is fractions and milliseconds and limits the resolution 
to select the minimum pause time during which the measuring current is applied to determine the UF at 
the p-n-junction. 

Recently, methods for determining the parameters of the electric-thermal model [15] are being 
developed, based on the analysis of frequency thermal characteristics, that is, the dependence 
of the amplitude and phase of change of the transition temperature TJ (œ) on the frequency when the device 
is heated with a power varying harmonically P(t)=P0+Pmsinot. In them, based on the calculation 
of the amplitudes and phases of the main harmonics of the heating power and temperature of the p—n 
junction, the module of transient thermal resistance (thermal impedance) and the @ phase shift between 
the temperature of the p—n junction and heating power is determined. Next, the dependence of the thermal 
impedance module on heating power modulation frequency Zthjc(@) and the thermal circuit parameters are 
determined. By analogy with electrical circuits, thermal impedance is considered as a complex number 
Zthjc=Zthjcejp, the real part of which determines Rthjc=Re[Zthjc ] thermal resistance. 

A technical solution: “Thermal Transient tester T3Ster” [16], which implements the methods of 
standards [11, 17] used to control the thermal resistance of semiconductor devices and integrated circuits is 
known. The device provides the ability to determine the direct voltage drop at the p-n-junction at 
a microsecond time resolution after the end of the heating pulse. Based on these data, using the software, 
the tested chip thermal resistance as a whole and of individual links of the thermal circuit are calculated. 
The TTR curve analysis obtained by any method is usually performed by classical methods, for example, 
by analyzing the first and second derivatives [12]. 

Recently, when analyzing the TTR curve, the apparatus of structure functions [16, 18-20] has been 
often used, which allows determining thermal resistance components. In accordance with it, at the first stage, 
using the inverse convolution of the graph of the first TTR derivative in the logarithmic time scale z=In(t) 
and the weight function wz=exp (z-exp(z)) the dependence of the spectral density of time constants (TCSD) 
is built. Next, it is discretized with a small dz step and a thermal model is constructed according to the Foster 
scheme. At the second stage, this scheme is transformed into the Cauer scheme, for example, according to 
the method in [2], and on its basis a set of dependences of the heat capacities on the thermal resistances of 
the elements included in the resulting scheme is obtained. This dependency is called the structure function, 
since for a one-dimensional heat flow, it characterizes the layers that make up the design of the sample under 
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study. This is manifested in abrupt changes in the graph of the structure function, and is easily determined on 
the graph of the first derivative of this function (the differential structure function). 

The method of structure functions is constantly refined and supplemented [21-24]. Critical factors 
for it are the noise level, which is always present in the TTR experimental determination and the need to 
analyze relatively small values of thermal characteristics, which may result in “peak emissions” in 
the beginning of the structural function graph, which in turn happens due to the influence of the first of 
the above mentioned factors. The use of well-known methods of filtering signals [16, 25] to reduce 
the influence of the first factor on one hand additionally requires the selection of a filter, and on the other 
hand leads to a distortion of the real picture of the occurring thermal processes. 

In [26-29], the authors used the original method of thermal relaxation spectrometry to analyze TTR. 
In [27], the thermal parameters of high-power bipolar transistors in plastic cases TO-252 and TO-126 were 
estimated using thermal relaxation differential spectrometry. For this purpose, a relaxation impedance 
spectrometer developed by the authors [28] was used, which makes it possible to determine the TTR. 
The thermal constants of the instruments and the distribution of the thermal resistance structure are 
determined in the form of discrete and continuous spectra. The continuous spectrumis calculated on the basis 
of higher-order derivatives of dynamic thermal impedance and corresponds to the Foster model, the discrete 
one corresponds to the Cauer model. The structure of the thermal resistance of the devices was represented in 
the form of a s1ix-lmnk electric-thermal RC model. 


3. METHODOLOGY AND MATERIALS 

For the experimental TTR determination, a method that does not require periodic connection of 
the test HFSVC IC [10] was tested. The crystal is heated by a direct current passing through a power diode 
located on the crystal and connected in parallel with the power transistor. Measurement of the crystal 
temperature is performed by measuring the forward UF voltage at another p—n junction at a small measuring 
current and then converting UF voltage to TJ temperature, using the TJ (UF). The resulting T(t) dependence 
can be approximated by several exponents. 


T, (t) = A,(1—e?"") +B, (1—e”"") (1 


The transient thermal characteristic can be determined using application software and a tester [9], 
which is designed to control the electrical parameters of high-frequency switching voltage converters 
(HFSVC) under mass production conditions. As a result of the UF(t) experiment, an array of digits is 
obtained, which should be converted into an array of digits T} (t) and then T(t), rae T(t)}=T>(t)-Ta, T is the 
heating temperature of the HFSVC crystal; T> 1s the total HFSVC crystal temperature; To is the ambient 
temperature. 

The next step is the process of approximation of the experimental TTR. Approximate TTR can be 
a different number of exponents. The authors used the Matlab tools [30]. To find the approximate TTR 
dependence in Matlab, the least squares method is implemented using the Leven-Berg-Marquardt 
optimization algorithm. To approximate it is necessary to perform the following operations: 

a. To convert the experimentally obtained matrix UF(t) to TJ(t), using a previously plotted temperature 
coefficient of voltage coefficient for a heat-sensitive parameter (TDP), which is one of the internal 
p-n-transitions of the microchip crystal. It also recalculates the heating power of the crystal temperature 
array to |W, taking into account the known value of the stabilized heating power. 

b. To set the type of approximating dependence. For example, for an up-proxy dependence for an 
independent time variable x based on six exponentials with 12 sought-for coefficients (A, B, C, D, E, F, 
pl, p2, p3, p4), its appearance is as follows: 


F(x) =A(l—-e *")+Bd-e ”)+Cd—e ’")+D0—-e"")+EFd-e ")+FU-e*) (2) 


To determine the lower and upper limit of the search range for each desired coefficient, which is 
included in the dependence according to step 2. To determine the initial approximation for each desired 
coefficient, which is included in the dependence according to step 2. To search for dependencies using 
the Matlab fit O function. To perform the analysis of the obtained dependency on additional output 
parameters of combinations of descriptive statistics produced by the fit function, such as: Determination 
coefficient R2 (rsquare) and matched R2 (adjrsquare), sum of squared errors (sse) and mean square error 
(rmse). Then complete the analysis with the help of visual inspection and validation of the graphs of 
experimental data and data according to approximating dependence by plotting these graphs in the same axis. 

The limits of the ranges in Section 3 can be selected from a range of positive integer values, 
for example, zero for the lower limit, and 1000 for the upper limit, or selected based on the values of 
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the coefficients obtained during the first iteration of finding the approximation. The initial approximation for 
each desired coefficient entering the approximating dependence can be chosen randomly in the range from 0 
to 1, or the values found during the first iteration of finding the approximating dependence. 


4. RESULTS 
For HFSVC IC LM2676 in TO-220-7 case in the absence of an additional radiator for heat removal, 
the time dependence ofthe crystal temperature is obtained to (3). 


T(t) =2,473(1— 6) $1, 125(1 2") +. 2, 773 (1 — 6 0) + 


(3) 
ais 24, 39(1 — gee: ) a 3, 366(1 _ Zs ) 4 23, 3(1 _ a ee ); 


Figure | shows the result of approximation by six exponentials. In the upper graph, the two 
dependencies experimental and approximating merged and are visually indistinguishable. The bottom graph 
shows the difference between these dependencies. 
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Figure 1. Graphical user interface of software toolkit in Matlab 


The obtained analytical dependence agrees well with the experimental one. This difference in 
the whole range does not exceed 0.26°C. According to the data of this equation, it is possible to construct 
a six-letter model in the form of the Foster scheme, in which the resistances are equal to the initial values of 
the exponentials, and the capacitances are determined by dividing the corresponding time constant by its 
resistance. So for the first RC element, we determine R1=2.473 Ohm, t=1/[p-1/0, 1794=5,574 c, 
C=1/R1=13.79 F. 

It is important to note that it is impossible to unambiguously compare each resistance to some kind 
of physical transition (for example, crystal-crystal holder (die), etc.), since the heat fluxes inside and outside 
the microcircuit are parallel, and the model is represented by a series connection of links. The standard [16] 
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and paper [14] also state that the FOSTER model ‘cannot be used to characterize the thermal structure. 
Real thermal capacitances are always connected to the base node-similarly to electrical “ground”, since 
the stored thermal energy is proportional to the temperature of one node, and not the temperature difference 
of two nodes, which is determined by the FOSTER model’[16]. This is a significant drawback of 
the FOSTER model, although it is a fairly correct mathematical model of transitional behavior. 

Similarly, approximations of four, three, two, and one exponent are obtained. These results are 
shown in Figure 2, in which the approximating dependence graph is shown in gray while the experimental 
data in black. Figure 2(a) shows the TTR approximation by four exponents. Figure 2(b) shows a version of 
the TTR approximation by three exponents. The obtained analytical dependences are in good agreement with 
the experimental ones. Mismatch graphs are visually indistinguishable. The temperature difference between 
the calculated and experimental TTR does not exceed 0.26°C. 

Figure 2(c) shows a variant of the TTR approximation by two exponents. The resulting analytical 
dependence T, (t) = 24, 64(1-e 8") 4.32, 5(1—e 8") has a mismatch in the initial part of the experimental 


TTR not more than 1,3°C. Figure 2(d) shows a variant of TTR approximation by one exponent. Visually 
dependencies mismatch. The resulting analytical dependence T(t) =56, 045(1—-e°"" ) has poor agreement 
with experimental results. The temperature difference between the calculated and experimental TTR 
in the initial part reaches 4.2°C. A significant deviation of the approximating dependence on the experimental 


one is observed in the initial part of the TTR. For this part, you can performan approximation again with two 
or three exhibitors as shown in Figure 3. 
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Figure 2. Approximation results by (a) Four exponents, (b) Three exponents, (c) Two exponents, 
(d) One exponent 
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Figure 3(a) shows a variant approximation by three exponents. The resulting analytical dependence 
T,(t) =2,35(1 —e") 40.5090 —e*?™ ) + 39,40. -—e °°") agrees well with the experimental results. 
The temperature difference between the calculated and experimental TTR does not exceed 0.05°C. 


This temperature difference is shown in the lower graphs of Figure 3. Figure 3(b) shows a variant of 
approximation by two exponents. The resulting analytical dependence T,(t)=10%1—e'?")+2890-e°"") 


agrees well with the experimental results. The temperature difference between the calculated 
and experimental TTR does not exceed 0.13°C. 
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Figure 3. Approximation results for the initial TTR par, (a) Variant approximation by three exponents, 
(b) Variant of approximation by two exponents 


It should be noted that the dependences for the initial TTR area have exponents that are significantly 
different from the total TTR exponent. One of the exponents has a time constant of less than one second, 
which characterizes the rapid change in temperature. On the basis of the obtained analytical dependencies, 
if necessary, it is possible to construct electrical models in the form of chains of RC elements in which 
the change in voltage U(t) during the transient process will coincide numerically with the TTR. Foster 
multilnk schemes are often used, in which each exponent corresponds to a link consisting of parallel- 
connected R and C. In this case, the initial value of the exponent determines the resistance R and the time 
constant t=1/|p| determines the capacitance C=t/R. 

Figure 4 shows the thermoelectric chip model, in which the current source J simulates the heating 
power P (thermal losses), the potentials at the nodes of circuit 1 and O simulate the temperature of the UJ p-n 
junction the ambient temperature Ta, the EMF source E imitates ambient temperature Ta. In this case, 
the total temperature of the p-n junction TJX(t)}=TatTJ(t) contains the constant component Ta 
and the variable TJ(t), due to heating. To simplify writing in the following equations, only the variable 
component of the temperature TJ(t) due to heating is considered. 





(a) (b) 


Figure 4. Representation of the electric-thermal model by, (a) The fosterdiagram, (b) The cauer 
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Obviously, one should not look for a correspondence between electrical and thermal resistances. 
Resistances R1, R2, R3 are calculated. Together with the capacitances C1, C2, C3 they provide a change in 
the voltage U (t) during the transient process (when the current source J is connected) numerically coinciding 
with the change in the crystal temperature (TTR) of the microcircuit. The resulting thermoelectric model can 
be implemented in any program for modeling electrical circuits (Orcad, Multisim, Matlab/Simulink, etc.) and 
investigate various transient conditions with regard to thermal processes. The voltage changes over time U1 
(t) will correspond to similar changes in the crystal temperature of the microcircuit T (t). 

In addition, the differences in the TTR circuits experimentally obtained in the heating mode of 
the sample and in the cooling mode, with the same thermal power, were also analyzed. These TTRs match 
quite well. The maximum discrepancies of the graphs were at a time level of 100s, on the order of 1—4 
for various suitable samples. For other time intervals of the same TTR there are even better coincidence of 
the graphs. 

Figure 5 shows as an example, an oscillogram of the passage of an electric pulse on a thermal 
electric model. Considering the oscillogram of the electrical transient process with a short-term change in 
the magnitude of the current source J, we can estimate the nature of the transition thermal process with 
a corresponding change in the heating power P. Thus, with the passage of a current pulse of 2J (double value) 
with a pulse duration of 60s, the voltage U1 changed from 57 V to 128 V, which corresponds to a change in 
the temperature of additional heating of the HFSVC crystal from 57°C to 128°C. At the same time, the total 
crystal temperature for an ambient temperature of 20 °C will change from 77°C to 148°C, which is a value 
close to the threshold value of the response of the temperature protection of the IMS VIP. A return to 
the nominal temperature of 57°C lasts about 1000 s. 

To verify the correctness of the model, a physical experiment was performed. The microcircuit was 
heated through a power diode for 5 minutes and on the free p—n junction the forward voltage UF was 
determined, which was later converted to the temperature TJ. The obtained dependence TJ (t) coincides well 
with the graph in Figure 5. 
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Figure 5. The passage ofheat pulse on the model 


5. DISCUSSION 

When forming the TTC electric-thermal model, the most important is the choice of the model 
structure, which in the problem under consideration corresponds to the determination of the number of 
exponents. In this case, it is often difficult to obtain an approximating dependency well coinciding with 
the TTC initial area. To reduce the influence of these factors when obtaining an approximating dependence, 
it is possible to use the spectraldensity of time constants (SDTC) [16, 19, 20, 23]. 

In [11, 20, 23, 25] it is proposed to do this by sampling a continuous graph of the spectral density of 
time constants with a fixed step dz. This, in turn, leads to an increase in the number of exponents and 
complication of the model, while simultaneously increasing the difficulty of calculation associated with 
the use of numbers with increased accuracy (double precision of numbers of the standard double float type 
with small dz is often not enough). 

The authors propose to use SDTC to determine the number of exponents, and to specify initial 
approximations when searching for constant time (coefficients at the exponent degree). The SDTC graphs of 
the experimental TTR, obtained by heating the LM2676 chip, is shown in Figure 6. The graph of the most 
important initial part of the structural function, obtained on the SDTC basis is shown in Figure 7. 

From the graph in Figure 6, four extremums are visible, which allows us to choose a model with 
four exponents. Additionally SDTC is a reliable source of information for specifying the initial 
approximation of the coefficients of exponent degrees, for which the choice of the initial approximation 
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and search ranges is most important for obtaining qualitative results in the approximation. As it can be seen 
from the SDTS graph, shown in Figure 6, it clearly identifies four time constants, the value of which 
practically coincides with the results of the approximation shown in Figure 2(a). 
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Figure 6. Graph of the spectral density Figure 7. Structure function graph 


of time constants 


Due to the initially inaccurate adjustment of the filters of the reverse convolution operation and 
the initial data of the experimentally determined TTR, the quality schedule of the SDTC is difficult to obtain 
at the first iteration. This is manifested in the fact that extremes can be shifted, which prevents you from 
immediately getting good approximation results (especially for the TTR initial part), or lowering 
the dimension of the problem, by reducing the number of desired coefficients in the approximation by 
eliminating constant time. The graph in Figure 6 was obtained by an iterative method of adjusting 
the deconvolution operation filter. 

In the method of approximation used in this work, the SDTC graph is used only to specify the initial 
approximations and ranges of the search for coefficients with exponential powers of the approximation 
dependence. This significantly improved the quality of the approximation while simultaneously reducing 
the number of exponents included in the approximating relationship because fromthe graph of the structural 
function as shown in Figure 7, the minimal set of these exponents does not explicitly follow and there always 
exists the probability of its inaccuracies due to filtering, as it has been mentioned already. This allows to 
obtaina minimized model of the thermal process in the IC without significant losses in the degree 
of its reliability. 


6. CONCLUSION 

The analysis of the methods for obtaining the thermoelectric IC model showed the main difficulties 
of their application for building models with a minimized number of RC circuits. The proposed solution is to 
refine the minimum structure of the thermoelectric model both on the basis of direct iteration based on 
the least squares method using the Levenerg-Marquardt optimization algorithm as well as using the SDTC 
graph as initial data, which allows reducing the influence of the filter factor in the deconvolution operation, 
with building a thermal model using the TTR structure function. The results obtained can be used to build IC 
behavioral models, taking into account the thermal processes occurring in them. 
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